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ABSTRACT 


Inverse FEM is a technique based on the loading force identification method of the Finite element. Charging 
approaches are based on the reverse finding problem approach. Linear, as well as non-linear system methods, are 
formulated based on the minimization of assumed objective functions. The least square error is primarily used as an 


analytical feature between the virtual and calculated device responses. 


The inverse finite element method is a hybrid numerical and experimental method to determine operating forces 
precisely. Mechanical members may not be designed considering the extreme cutting forces. Hence the problem is taken 
to identify the cutting forces of a cutting lathe tool during turning operation using inverse FEM technique. Strain gauges 


are mounted on a lathe tool near the fixed end locations to measure the strains during turning operation. 


Lab view are called virtual instruments is used to evaluate the impact forces on a Lathe tool. The Lab view 
includes a wide range of tools to evaluate, monitor and store data and tools to help you repair the coding. The turning is 
carried out with the instrument and strains are measured at different speeds at known locations. The determined stresses 


of a device are compared by a least square approach with the empirical strains. 


Studies have been carried out to determine impact forces using a Lathe tool employing Inverse FEM. This 
method is subsequently employed in the determination of cutting force during turning operation. Inverse FEM with strain 
gauge as sensors is used for estimating the cutting force for different speeds and for a speed of 250 rpm the cutting force 
estimated was 225.63N. The same force has been measured using a lathe tool dynamometer and measurement of force 
has found to an error of 15%. It is found that for higher speeds the error in measurement is reduced. i.e. at higher speeds, 


the vibrations are less and steady state is reached. 


The cutting force estimated by inverse Fem using strain gauge as sensors is compared with the force measured 
using lathe tool dynamometer. The error in measurement is found to be 15%. This is observed to reduce as the speed 
increases. During this study Inverse finite element approach is applied for estimating the impact forces on cantilever 


beam and operating forces in a lathe. 
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INTRODUCTION 


For the performance of certain tasks, structural members and machines are designed and developed. They are 


especially constructed on the basis of loads and boundary conditions taking into account the material's properties. 
Static load 


Static loading Static load does not change with time. These are loads that progressively increase over time or have 
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insignificant dynamic effects. Since static loads are much easier for structural analysis than dynamic ones, construction 


codes usually require statically-equivalent loads for wind, traffic or earthquake dynamic loads. 
Dynamic load 


The dynamic load changes over time. If this changes gradually, the response of the system can be calculated by static 
analysis but if the response changes quickly, then the dynamic analysis will decide the response. Because the system 
cannot respond to loading rapidly, a dynamic load can have a significantly greater impact than a static load of the same 


magnitude. Significant dynamic effects show dynamic loads. E.g. impact loads tides, wind gusts and extreme earthquakes. 
Inverse FEM 


Inverse FEM is a technique based on the loading force identification method of the Finite element. Charging approaches 
are based on the reverse finding problem approach. A number of approaches have been developed for the linear system in 
this field. Linear, as well as non-linear system methods, are formulated based on the minimization of assumed objective 
functions. The least square error is primarily used as an analytical feature between the virtual and calculated device 
responses. The identification of model parameters is crucial for different algorithms for resolving reverse problems such as 


the identification of Force. Rapid numerical device creation and computer technology growth make Inverse FEM possible. 


General Procedure of Solving an Inverse Problem 


Inverse analysis 


General inversions, optimization 


Regularization 


Figurel: Flow Chart to Solve Inverse problems 


Impact Factor (JCC): 10.2746 NAAS Rating: 3.11 


Estimating the Impact Forces and Operating Forces in a Lathe Tool Using Inverse Finite Element Approach 81 


STATEMENT OF PROBLEM 
Application Problem on Lathe Tool 


68.863 Modi fed parameters 
= | FEM System 
ss | Resuts of 
- bad > Simulation 
Ah 9 . 
2 se c wi 
» 2° ro iser ; Compare with 
oe 55 Module meesured , 
parameters Ow Por 
Front view ——— . 
Seale: 1:1 Inverse FEM model 
Pocerct ect Ome 


Mechanical members may not be designed considering the extreme cutting forces. Hence the problem is taken to identify 
the cutting forces during turning operation using the inverse FEM technique. As we know lathe tool acts as a cantilever 


beam. The inverse identification problem on cutting lathe tools is carried out in the laboratory. 


A model of a cutting lathe tool from industrial application is considered for our laboratory experimental work. 


Strain gauges are mounted on a lathe tool near the fixed end locations to measure the strains during turning operation. 


MATHEMATICAL THEORY ON INVERSE PROBLEM 
Formulation and Solution of Inverse Problems 
General Inversion of System Matrix 
Consider a system whose forward problem can be formulated as Equation 3.2. For the calculation of input X, the 


solution to the related inverse question can be sought. 
e _ = 2 m 
X yx1 = SyxuYuxi 


Where the superscripts ' e' and' m' are respectively the value calculated and determined. The output vector Y is 
now known as the device results obtained by means of regular measurements, documents or observations, and S is an 


opposite of matrix S generated by using Equation 3.31. 
¢ Under-posed problems- when M<N 
¢  Even-posed problems- when M=N 
* Over-posed problems- when M>N 


The following section describes the way in which these reverse problems are solved. While assuming that S can 
be obtained and thus X can be calculated by using Equation 3.31. The quality of the prediction following the reversal must 


be checked. Production Y can be calculated with approximate X, based on the forward Equation 3.2 formula. 
Y = SuxwX wer = Is mw S xu Mua (2) 
This equation clearly shows that, if 
IS wah | =f (3) 
Where I am a matrix marker, calculation data can be repeated by the reverse process. When a S-g obtained does 
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not satisfy Eq 3.33, it will not replicate the measurement results. The matrix 


Roux) a Is MXN NEM (4) 


Therefore, defined as the matrix of reproductivity output. In terms of reproducing input or output data, the 
computer Ro may give an indication for the accuracy of an inverse operation. If Ro= I, the reversal is replicable. Next, the 


real X is claimed to fulfill the forward formula for the measured Y obtained for the actual event: 


Ve 


MX1 — Si (5) 


Where the superscript t stands for the true values. Substituting this equation into equation 3.31 gives, 
e = t 
X nxt = |S585) Sane ee (6) 


The preceding equation clearly shows that, if 


een | =I (7) 


The inverse procedure can provide the exact estimation of the inputs. Therefore, the matrix 


Riyyw) = SiS i | (8) 
Is defined as an input reproducibility matrix. Computing Ri can give an indication of the quality of an inverse 
procedure in terms of estimating the input of the system. 
Under-Posed Problems: Minimum Length Solution 


If M is less than the estimated number of unknowns for the system, then it is a fundamental problem. There is an_ infinite 
number of solutions to the underlying problem that follow the exact equation without the defect of the accompanying 


forward model. 
Even-Posed Problems: Standard Inversion of Matrix 


If the number of known results is equal to N, the problem is even and the device matrix SN x N square matrix with rows an 


d columns of N is the same. The unknown number is calculated for the program. The machine has the same number. 
Over-Posed Problems: Least-Squares Solution 


If M, which is larger than the unknown unit, N, the calculated number of unknowns is over imposed. There can be a 
number of solutions to the problem of over position to satisfy certain right forward model equations. A solution which is 
physically significant and properly compromises all the equations of the future model, therefore, needs to be found. The 
"Lesser Squares" approach is a common method. The LS solution is formulated in the following way. DefineITas a first 


function. 
II =(Y™— SX)™(Y™-SX) (9) 


In addition, this is the L2 rule for the estimation and measurements of the forward model. 
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oll _y (10) 
Ox 
oll __ssr(ym_sx)=0 (11) 
Ox 
Or 
STy™ = S™SX (12) 


Consider the matrix S'S is symmetric as 

(STS) T +(S) 7 (S") T= STS (13) 
The approximate X can be calculated as if STS were invertible. 

x =(S's/' Sty" (14) 
The definition of the generalized inverse matrix for the over-posed inverse problem is 

S=[STS]" st (15) 
The output reproducibility matrix is 

Ro=SS%=S([S'S]?S? (16) 


Usually, this cannot be a matrix of identity. Consequently, generally speaking, the LS solution will 


be reproducible. The reproducibility matrix of the input is, however, 
ReS*S=(S'S]'S'*S=I (17) 
This implies that the LS solution is input reproducible and gives the true estimation of the input of the system. 


Lathe Tool model 
Measurement of Cutting Force on a Lathe Tool 
The turning is carried out with the instrument and strains are measured at different speeds at known locations. The 


determined stresses of a device are compared by a less square approach with the empirical strains. 
Forward Analysis 


The tool Lathe is Catia V5 modeled. For meshing the model, the element Solid45 is used. The meshed tool is simulated by 
limiting the whole D.O.F to one end of the tool in ANSYS as a cantilever strap and applying the statically loading range to 


the open end of the tool. Strain is documented at known places. 
The material and geometrical properties are, 
E=2.1x10°N/mm/?; u=0.3 


A =238.177mm’, p=7850e-09 N/mm? 
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Figure 4: FE Model with Boundary Conditions Figure 5: Strain Plot for Applied Force 


ANSYS Strain Readings for Range of Static Load on Lathe Tool 


Static Load vs strain curve 
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Figure 6: Static Loads vs. Strain Curve 


EXPERIMENTAL MODEL 


A cutting lathe tool is taken for experimentation in the laboratory and Strain gauges are mounted at predetermined points 


on the surface of the lathe tool. The cutting tool is operated on a lathe for different speeds by maintaining the constant 
depth of cut (1.5mm) and constant feed (0.12rev/min). 


Figure 7: Location of Strain Gauge on Lathe Tool Figure 8: Experimental Set Up 


Table 1: Experimental Strain Readings for Lathe Tool Operated at 
Different Speeds 


250 15.—C«dLSC (iC 15.47 — 
350 15 0.12 108.5 
500 15 0.12 168.5 
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strains (microns) 


Lathe tool operated at 250rpm speed 


_——strains at 250 rpm speed 


time (milli seconds) 


Lathe tool operated at 350 rpm speed 


——strains at 350 rpm. 


strains (microns) 


time (milli seconds) 


Experimentally Measured Strain Readings at 250 Rpm | Experimentally Measured Strain Readings at 350 Rpm 


strtains (microns) 


Lathe tool operated at SOOrpm speed 


strains at S00 rpm speed 


time (milli seconds) 


Experimentally Measured Strain Readings at 500 Rpm 


Table 2:Least Square Error Readings for Cutting Speed 250 Rpm 


206.01 68.5 75.47 485 
215.82 71.8 75.47 134.68 
225.63 75.1 75.47 1.369 
235.44 78.3 75.47 80.0 
245.25 81.6 75.47 375 


Table 3: Least Square Error Readings for Cutting Speed 500 Rpm 


490.5 163.21 168.5 27.9 
500.31 166.47 168.5 4.12 
510.12 169.73 168.5 1.51 
519.93 173.00 168.5 20.25 
529.74 176.26 168.5 60.21 
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Least square plot for cutting speed 250 rpm least square plot for cutting speed 500 rpm 


200 —least square plot —least square plot 


Minimum error(*e-13) 
w 
S$ 
So 

Minimum error (*e-12) 


load (N) Load (N) 


Plot of Least Square Difference for Cutting Speed 250 Plot of Least Square Difference for Cutting Speed 
Rpm 500 Rpm 


Table 4: Cutting Force Readings from Lathe Tool Dynamometer 


250 1.5 0.12 264.87 
350 1.5 0.12 372.78 
500 i Ges) 0.12 549.36 


Table 5: Comparison of Cutting Force from Inverse FEM and Tool Dynamometer 


350 323.73 372.78 13.15 
500 510.12 549.36 7.14 
RESULTS AND DISCUSSIONS 


Studies have been carried out to determine impact forces using a Lathe tool employing Inverse FEM. This method is 
subsequently employed in the determination of cutting force during turning operation. Following is the summary of the 


results obtained during this study. 
Lathe Tool 


The methodology developed above is applied for measuring force on a lathe tool during turning operation. Inverse FEM 
with strain gauge as sensors is used for estimating the cutting force for different speeds. Strain gauges are mounted on the 
cutting tool so as to measure the cutting forces in the vertical direction, and for a speed of 250 rpm the cutting force 
estimated was 225.63N. The same force has been measured using a lathe tool dynamometer and measurement of force has 


found to an error of 15%. 


The experiment is repeated for different cutting speeds and it is found that for higher speeds the error in 


measurement is reduced. 1.e. at higher speeds, the vibrations are less and steady state is reached. 
CONCLUSION AND SCOPE FOR FUTURE WORK 


For optimal design of the machine and structural members, it is necessary that the forces to which such members are 
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subjected to be known before the design process is started. But usually, the operating forces are not exactly known and 


design based on such approximate loads is not optimal. This shortcoming is overcome by employing the inverse finite 


element method. The inverse finite element method is a hybrid numerical and experimental method to determine operating 


forces precisely. 


The method is extended to measure the cutting force in the lathe tool. The cutting force estimated by inverse Fem 


using strain gauge as sensors is compared with the force measured using lathe tool dynamometer. The error in 


measurement is found to be 15%. This is observed to reduce as the speed increases. 


Thus, during this study Inverse finite element approach is applied for estimating the impact forces on cantilever 


beam and operating forces in a lathe. 


Future Work 
1. The work can also be extended to find the dynamic forces acting on the automobiles and machine parts. 
2. The inverse method can also be carried out to find the Material property and Geometry identification. 
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